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ABSTRACT 

Accurate photometric redshifts are calculated for nearly 200,000 galaxies to a 4.5 micron flux limit of ~ 13 
/xJy in the 8.5 deg 2 SpitzerflRAC Shallow survey. Using a hybrid photometric redshift algorithm incorporating 
both neural-net and template-fitting techniques, calibrated with over 15,000 spectroscopic redshifts, a redshift 
accuracy of a = 0.06(1 +z) is achieved for 95% of galaxies at < z < 1.5. The accuracy is a = 0.12(1 + z) 
for 95% of AGN at < z < 3. Redshift probability functions, central to several ongoing studies of the galaxy 
population, are computed for the full sample. We demonstrate that these functions accurately represent the 
true redshift probability density, allowing the calculation of valid confidence intervals for all objects. These 
probability functions have already been used to successfully identify a population of 5/?/fzer-selected high 
redshift (z > 1) galaxy clusters. We present one such spectroscopically confirmed cluster at (z) = 1.24, ISCS 
J1434. 5+3427. Finally, we present a measurement of the 4.5/im-selected galaxy redshift distribution. 

Subject headings: galaxies: distances and redshifts — galaxies: clusters: individual (ISCS J 1434.5+3427) 
— galaxies: clusters: general — galaxies: evolution — galaxies: photometry — methods: 
statistical 



1. INTRODUCTION 

In modern wide-field imaging surveys, accurate photomet- 
ric redshifts have become an indispensable tool for studying 
the evolving galaxy population. This approach allows stud- 
ies of the galaxy luminosity and correlation functions, power 
spectrum analyses, and rare-object searches, reserving costly 
8-m class spectroscopic follow-up for analyses where such 
data are crucial. 

Much of the new extragalactic science enabled by the 
Spitzer Space Telescope will not just be more efficient with 
photometric redshifts, but rather completely reliant on them. 
Spitzer is uniquely capable of detecting huge numbers of 
galaxies from their rest 1.6 fim stellar emission or mid-IR 
PAH features all the way to z = 2 and beyond. These objects 
are often quite faint in the optical, due to either quiescence or 
heavy extinction, making spectroscopic follow-up difficult or 
impossible with existing faciliti es. 

The IRAC Shallow survey fcisenhardt et alJ 12004 is a 
wide-field 8.5 deg 2 Spitzer/IRAC imagin g survey in the 
NOA O Deep Wide-Field Survey (NDWFS; Jannuz i and Devi 
119991) Bootes field, designed to study galaxy formation and 
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evolution across a wide range of redshifts, mass scales, 
colors and environments. Extensive complementary data 
includ es a 5 ksec per pointing X-ray mosaic (XB ootes; 
iMurrav et alJ2005tlKenter et aljioollBrand et alJ2006l) . and 
SpitzerMWS 24/70/160 /zm imaging (GTO program, Soifer 
PI) over the full survey area. In addition, the FLAM INGOS 
Extragalactic Survey (FLAMEX; Elston et al. 2006) provides 
deep J and K s imaging over the central « 4.7 deg 2 . Finally, 
there are 17,017 spectroscopic redshifts of galaxies to z ~ 0.8 
and AGN to z = 5.85 from the A GN and Galaxy Evolution 
Survey (AGES: ICool et al J 120061 Kochanek et al. in prepa- 
ration), along with ~ 500 redshifts of fainter, higher red- 
shift galaxies from various other surveys in Bootes, which are 
available for calibration of the photometric redshift algorithm. 

In this paper, we combine the optical photometry from the 
NDWFS, the near-IR photometry from FLAMEX, and the 
mid-IR photometry from the IRAC Shallow survey to com- 
pute accurate photometric redshifts for 4.5 /im-selected galax- 
ies to z = 1 .5 and AGN to z = 3. A hybrid method, combining 
a standard template-fitting algorithm with an artificial neural 
net-based approach, is developed to make optimal use of the 
large AGES spectroscopic calibration sample. 

Redshift probability functions are derived for the full sam- 
ple, nearly 200,000 objects over 8.5 deg 2 , and are employed 
to calculate the 4.5/xm redshift distribution. They are also the 
input to a high-redshift cluster search underway in Bootes, 
which has already y ielded the highest r edshift galaxy cluster 
to date ((z) = 1.41; Stanford et al. 2005). This paper presents 
a new, spectroscopically confirmed high redshift cluster at 
(z) = 1.24. 

This paper is organized as follows. In ^2] the various data 
sets used in this paper are briefly described. In <j3]the photo- 
metric redshift algorithm is presented, and the results are com- 
pared to the AGES spectroscopy. The calculation of redshift 
probability functions is described in and in |5] we employ 
them to compute the 4.5/xm redshift distribution and present 
a new (z) = 1.24 galaxy cluster discovered using them. We 
summarize our results in |6] All magnitudes are Vega-based. 
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2. PHOTOMETRIC AND SPECTROSCOPIC DATA 
2.1. Optical and NIR Photometry 

Optical B w , R, and / imaging data were taken from the 
third data release (DR3) of the public NDWFS survey in the 
Bootes field. These data, obtained with the Mosaic-I cam- 
era on the NOAO 4-m, are described fully in Jannuzi et al. 
(in preparation) and available through the NOAO Science 
Archive ( http://archive.noao.edu/nsa/). 

Robust photometric errors were estimated via Monte-Carlo 
simulation, and extensive flagging of both pixel and image 
artifacts allows selection of high-quality photometric samples 
over the full survey area. The optical photometry reaches 3a 
5" diameter Vega depths of B w = 25.3, R = 24. 1, and / = 23.6. 
The large apertures are taken to better match the IRAC PSF 
and to minimize the effect of seeing variations in the optical 
data, as discussed below. The NDWFS is significantly deeper 
for point source or small aperture photometry. 

The FLAMEX survey JElston et alJ l2006) is a near-IR J- 
and Tfj-band imaging survey undertaken with the FLAMIN- 
GOS camera on the NOAO 2.1-m. Photometric errors were 
determined via extensive Monte Carlo simulations, account- 
ing for PSF variations across the field. The survey covers the 
central ^4.7 deg 2 subset of the Bootes field to a 5" diameter 
aperture, 50% completeness limit of K s = 19.4. 

2.2. IRAC imaging 

The IRAC Shallow survey, introduced in Eisenhar dt et alJ 
(2004), is a Spitze r/IRAC imaging survey in the NDWFS 
Bootes field, covering 8.5 deg 2 with 3 or more 30 second ex- 
posures per position at 3.6, 4.5, 5.8, and 8.0 /im to 5a flux 
limits in a 5" aperture of 10.0, 13.3, 78.0, and 68.3 ply, or 
Vega magnitude limits of 18.6, 17.8, 15.4, and 14.9, respec- 
tively. 

Separate photometric catalo gues were extracted in ea ch 
channel using SExtractor 2.3.2 (Bertin and Arnouts 1996) in 
double-image mode, producing matched catalogs in the other 
3 IRAC bands. The detection images in each channel were 
weighted by the e rror images generated by the MOPEX mo- 
saicking software (Mako vozand Khadll 2005). This paper fo- 
cusses on the 4.5 ^m-selected catalog, which is the natural 
selection band for the 1 < z < 2 cluster search described be- 
low (e.g., Eisenhardt et al., in preparation). 

Quality control was maintained through the extensive use 
of flags. In particular, to ensure reliable colors, objects with 
less than 3 exposures in any single pixel within the aperture of 
interest were rejected from the final catalogs. This primarily 
removes objects from the edges of the field, with only ~ 1 % of 
the 4.5 pm sample rejected in the main overlap region. This 
spatial selection function is well quantified and will not af- 
fect any science analyses. This flagging was carried out for 
all of our apertures, which span diameters of 1-20", chosen 
to match the DR3 NDWFS catalogs that form our primary 
complementary data set. All IRAC aperture photometry was 
corrected to large (24.4" diameter) apertures to account for 
PSF losses. For the 5a 5" diameter aperture of interest in this 
work, there are 21 1, 260 objects which have the full exposure 
time in both the [3.6] and [4.5] bands. 

2.2.1. Catalog Matching 

There is a small offset, of 0.38" in right ascension and 0. 15" 
in declination, between the astrometric solutions of the near- 
and mid-IR catalogs, which are tied to the 2MASS reference 
frame, and the optical catalog, which is tied to the USNO- 



A2 (M onet et alJl200l . This offset, caused in part by errors 
in centroiding the bright Tycho-II stars used to zeropoint the 
USNO-A2 astrometry, was removed prior to matching (for 
further details, see Jannuzi et al., in preparation). 

Detections in the optical and near-IR were matched to the 
4.5 /im sources if the centroids were within 1" of each other. 
For extended objects, detections in the different bands were 
matched if the centroids were within an ellipse defined us- 
ing the second order moments of the light distribution of the 
object 12 . IRAC-selected objects with no match in a given op- 
tical or near-IR band were assigned a Monte Carlo estimated 
la flux limit representing the sky variation in a 5" diameter 
aperture. 

2.2.2. Photometric Redshift Sample 

We define the photometric redshift sample as the subset of 
the [4.5]-selected matched catalog for which both the [3.6] 
and [4.5] data have the full 90s exposure time, and for which 
at least 2 of the 3 optical bands contain useful (i.e., unmasked) 
photometric data. Note that object detections in the optical, 
near-IR, and [3.6] bands are not required; la limits are used 
for non-detections. This results in a final multi-wavelength 
photometric redshift sample of 194,466 objects. 

2.3. Spectroscopic Redshifts 

The large sample of spectroscopic redshifts from AGES 
(Kochanek et al. in preparation) provides a crucial train- 
ing sample for the photometric redshift algorithm. AGES 
is a wide-field MMT/Hectospec (Fabr icant et al.ll2005l) red- 
shift survey, version 2.0 of which contains high quality spec- 
troscopic redshifts for 17,017 objects, including galaxies to 
z ~ 0.8 and AGN to z = 5.85 dCool et alJl2006l) . From this 
sample, 15,052 objects correspond to sources in the photo- 
metric redshift sample defined above. 

The AGES survey was designed to allow magnitude- 
limited samples to be selected as a function of wavelength 
from the X-ray to the radio. In the IRAC 4.5pm band, the 
AGES survey is w 100% complete to 15.2 mag (150 /iJy), and 
statistically complete (with 30% random sampling) to 15.7 
mag (95 pJy). 

However, more than half the spectroscopic sample lies be- 
yond the fainter of these limits, driven primarily by the optical 
magnitude limits of / = 20 for extended objects and / = 21.5 
for point sources, where the latter limit was designed to in- 
clude relatively large numbers of AGN and quasars. In addi- 
tion, the complete samples in the X-ray, UV, Far-IR and radio 
result in the survey being overweight in both active galaxies 
and those which are strongly starforming. 

We have also assembled a deep (R ~ 25) heterogeneous 
sample of spectroscopic redshifts from several ongoing 
projects in the Bootes field for the purposes of calibrating the 
photometric redshifts. These projects include spectroscopic 
studies of optical, near-IR, IRAC and MlPS-selected sources, 
including our own spectroscopy of 1 . 1 < z < 1 .4 galaxy clus- 
ter members described below. As many of the principal in- 
vestigators of these projects are also co-investigators of the 
present work (AD, BTJ, BTS, DS, SAS), we will call this the 
"in-house" sample for brevity. 



12 This ellipse was defined with the SExtractor parameters 2 X 
A_WORLD, 2 X B_WORLD, and THETA_WORLD. 
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3. METHOD: PHOTOMETRIC REDSHIFTS 

The impressive accura cy of photometric redshift algo- 
rithms of various typ es (Bemtez 2000; Brodwin et alJl2006t 
Connolly etallll997t iFernandez-Soto et alJl2002t iFirth et al] 
2002J Fontana et al] 120001 iHsieh et al.1 l2005t ISawicki eTaLl 
19971 IWolf et al.N2003l) is a good indication that the field is 
rapidly maturing at least as concerns optical/NIR surveys. It 
is less clear which methodology is best for redshift estima- 
tion from photometry extending into the L and M bands and 
beyond. A host of physics not included in current popula- 
tion synthesis models, including ubiquitous PAH emission 
and molecular and silicate absorption, could potentially com- 
plicate photometric redshift estimation in the first generation 
of Spitzer surveys. 

Well-calibrated empiri cal me t hods such as f itting func- 
tions (e.g. iBrunner etal] 120001 IConnoll vet all 119951) and 
neural nets jJCollister and Lahavl 12004 IFirth et alJ boolt 
Vanzella et al. 2004) are capable of matching or surpassing 
the accuracy of template-fitting methods. They can be par- 
ticularly powerful for analyses in which it is sufficient to pre- 
dict accurate redshifts for a relatively small subset of the gen- 
eral galaxy population, albeit one with a l arge, represe ntative 
spe ctroscopic training set. For instance. IB lake etal] (2006) 
and Padmanabhan et al. (2006) measured the galaxy power 
spectrum using only bright, low-redshift, early type galaxies, 
selected with a variety of color and magnitude cuts. The large 
4000 A breaks and lack of SED ambiguity produces highly 
accurate redshifts in these samples. 

In this paper we wish to develop a methodology with opti- 
mal photometric redshift accuracy over a much larger redshift 
baseline and spanning all galaxy types. Not surprisingly there 
is no spectroscopic sample spanning the full range of galaxy 
color, magnitude, spectral type, and redshift in the present sur- 
vey. Despite the large size of the AGES spectroscopic sample, 
it consists of the low redshift, high luminosity component of 
the general IRAC photometric redshift sample. The situation 
for QSOs and AGN is more encouraging. The broad range of 
targeting methods and deeper limiting magnitudes likely lead 
to a broadly representative sample of unobscured AGN for all 
redshifts. 

On the other hand, template-fitting methods calibrated with 
even modest spectroscopic redshift samples produce photo- 
metric redshifts which are generally robust outside of the nar- 
row parameter space in which they are specifically validated. 
For this reason template-fitting algorithms are the de facto 
standard in the literature. An important advantage of this 
technique is the straightforward generation of redshift prob- 
ability functions from the likelihood analysis, which are key 
to many science applications. This method is only accurate if 
the galaxy templates are representative of the observed galax- 
ies. Strong PAH-emitting galaxies, for which reliable tem- 
plates do not yet exist, along with quasars and AGN which 
have minimal continuum breaks, present a challenge for this 
technique. 

In the IRAC Shallow survey we therefore employ a hybrid 
approach in which a template-fitting algorithm is used as the 
core method, supplemented by a well-calibrated neural net 
technique for the small subset of objects which both require 
and merit it. 

3.1. Empirical Template-Fitting Algorithm 

The template-fit ting algorithm closel y follows that de- 
scribed in detail in iBrodwin et all J2006T) . Coleman, Wu & 



Weedman (1980, CWW) galaxy SEDs, supplemented by the 
iKinnev et all ill 9961) empirical starburst (SB3 and SB2) SEDs 
are used as basis templates. Th ese templates were ex tended 
to the far-UV and mid-IR using B ruzual and Charlol ( 120031) 
models. To improve the redshift accuracy linear interpolates 
were derived, resulting in 19 templates finely spanning the 
template space between the CWW Elliptical and SB2. 

Population synthesis codes do not yet accurately model 
the complicated physical processes, in particular strong PAH 
emission, which can dominate the rest-frame A > 5 /im emis- 
sion from normal galaxies. In view of this, and given that 
only the brightest ~ 10% of our sample has well-measured 
fluxes in the 5.8 and 8.0 /im bands, we elected to limit the 
fitting to the A < 5 /im regime where the models are expected 
to better approximate the true SEDs. This mild restriction still 
permits us to sample the stellar peak at 1.6 Mm, a u s eful red- 
shift indicator (e.g., ISimpson and Eisenhardtlfl999t ISawickil 
120021) out to z ~ 2. The PAH emission at A ~ 3.3 /im, which 
could potentially cause a template mismatch for strongly star- 
forming galaxies, is quite modest in terms of equivalent width 
(EW = 0.02um: lLu et alJ2003l) . Indeed, this feature has only 
0.5% of the powe r of the PAH features longward of 5/im 
llHelou et all200Cl) . Since it redshifts out of the [4.5] band by 
z ~ 0.5, any deleterious effect on redshift estimation should 
be limited to low redshifts. 

Redshifts were fitted between < z < 5 using 5" diame- 
ter aperture BwRI+ JK S + [3. 6][4.5] photometry where avail- 
able. Reliable detections are not required as this would im- 
pose a strong selection effect; limits are used for areas of sky 
which were observed but for which no object was detected. 
Many deeper IRAC surveys suffer from considerable confu- 
sion in the bluer bands, leading to substantial difficulties in 
deriving accurate photometry. Due to the combination of area 
and depth targeted in the Shallow survey, the images are not 
confused, allowing straightforward aperture photometry. 

To obtain robust galaxy colors, it is common in ground- 
based imaging surveys to smooth the images in all bands to 
a common worst seeing. Due to the relatively large PSF mis- 
match between Spitzer (~2") and the ground-based optical 
data (0.8"- 1.3"), it is not clear that this is the best approach. 
The large apertures required to enclose a substantial fraction 
of infrared light, empirically 5", are sufficient to minimize 
the effects of seeing variations in the optical and near-IR im- 
ages. We have verified that photometric redshifts computed 
with smaller apertures (3") exhibit redshift- and position- 
dependent systematic errors, likely due to ground-based see- 
ing variations. These effects vanish with the larger 5" aperture 
photometry. 

3.2. Correcting Templates and Zero Points 

Comparison of preliminary photometric redshifts with 
AGES spectroscopy indicated an error in the mid-IR color 
of model elliptical galaxies. The spectral synthesis codes at 
A > 1 /im model a zero-color Rayleigh-Jeans tail typical of 
simple st ellar po pulations. Recent Spitzer studies of nearby 
galaxies (Pahre et al. 2004) indicate that, while reasonable for 
late type galaxies, early-type galaxies are in fact blue in the 
mid-IR, with [3.6] -[4.5] « -0.15, due to CO absorption. 
To account for this the elliptical template (in Fx) was scaled 
down between 1 < A(/im) < 5 by a factor of m{\{pm) - 1), 
where the slope, m = -0.11, was determined by maximizing 
the photometric redshift accuracy for high S/N AGES ellip- 
tical galaxies. This change was carried through the template 
interpolations. 
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The large AGES spectroscopic sample was also used exten- 
sively to analyze the absolute inter-survey photometric cal- 
ibration. For those objects well characterized by the model 
templates, the secure spectroscopic redshifts allow the AGES 
sample to be effectively used as spectrophotometric standards 
to determine inter-survey photometric offsets. These were 
found to be negligible in the optical and small in the near- 
IR, AJ = -0.10 mag and AK S = -0.02 mag. However, in the 
mid-IR the offsets between the observed photometry and the 
iBruzual and Chariot! (120031) models were considerably larger, 
A[3.6] = 0.27 mag and A[4.5] = 0.32 mag. 

While systematic zeropoint errors at the ~ 5 — 10% level 
are possible with IRAC data jReach et al.ll2005l) . this much 
larger error is likely due to the inadequacy of the spectral syn- 
thesis mod els in the near- an d mid-IR. Indeed, as illustrated 
recently by Maraston (2005), the near- and mid-IR to opti- 
cal c olors pr edicted by independent spectral synthesis models 
(e.g., Bruzual and Chariot 2003; Fioc and Rocca-Volmerange 
119991 IVazdekis et all 1 19961 I Vazquez and Leithererl 120051) 
have a scatter of 0.2-0.3 mag even for identical input stel- 
lar evolutionary t racks. In addition , the new population syn- 
thesis models by Maraston (2005), which include the con- 
tributions from the post main sequence evolutionary phases, 
predict higher infrared fluxes than previous models, which 
are more consistent with our corrected aperture photometry. 
Work is in progress (Kochanek et al. in preparation) to em- 
pirically derive low resolution rest-frame mid-IR SEDs for 
both galaxies and AGN using the AGES spectroscopic sample 
and the full multi-wavelength photometry in the Bootes field. 
For the present paper we adopt the above magnitude offsets 
as empirically motivated corrections to the photometry-SED 
combination that optimize photometric redshift accuracy. 

This calibrated template-fitting algorithm was tested on 
AGES objects with z < 1, as all z > 1 AGES objects are 
QSOs or AGN and of little use in evaluating the performance 
of galaxy template-fitting algorithms. As shown in the top 
panel of Figure the results are quite poor, with an rms dis- 
persion in the photometric redshifts about the true redshifts 
of a « 0.38, or a « (1 +z)0.27. The dispersion for a 95% 
clipped sample, the 95% of objects with the smallest absolute 
redshift difference, is significantly better, with a w 0.09, or 
a « 0.07 (1 +z). While indicative that in the mean the method 
works quite well, there are clearly objects not well fit by the 
empirical templates. As discussed above, PAH-emitting and 
active galaxies, two object classes not represented in the em- 
pirical galaxy templates, are over-represented in the AGES 
sample. To see how this affects the statistics, the lower panel 
of FigureHomits objects likely to be either PAH-emitters or 
AGN. The former are defined in this work as objects with an 
IRAC color of [5.8]- [8.0] > 1. The latter were identified 
on the basis of their infrared colors in an IRAC color-color 
diagra m similar to the one described in detail in IStern et all 
ll2005l) . That paper demonstrated that active galaxies can be 
rea dily identified bas ed solely on their infrared colors (also 
see lLacv et aTl 120041) . Omitting these two populations re- 
moves 70% of the variance, or about half the rms error, while 
maintaining similar clipped statistics. We now turn to tech- 
niques of estimating redshifts for these populations. 

3.3. Artificial Neural Net Algorithm 

Although the utility of purely empirical methods is lim- 
ited to the parameter space of the calibration sample, such 
methods do offer several unique advantages. Neither abso- 
lute band-to-band calibrations, nor a complete knowledge of 
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FIG. 1. — Comparison of calibrated template-fitting photometric redshifts 
with AGES spectroscopy. In the top panel all z < 1.1 AGES objects are 
plotted, whereas in the bottom panel likely AGN and strongly PAH-emitting 
galaxies have been removed. The full rms dispersion markedly improves in 
the bottom panel. The clipped statistics (in parentheses) are quite good in 
both cases, with a < 0.1. 



the galaxy population are required for accurate redshift esti- 
mation. Given a large spectroscopic training set, algorithms 
such as polynomial fits or artificial neural nets (ANNs) can be 
trained to predict the redshifts of objects of all types using the 
actual survey photometry. 

Regardless of method, the uniqueness of the mapping from 
colors to redshifts is the underlying limitation in redshift ac- 
curacy. Objects with approximately power-law spectra such 
as QSOs will never allow very accurate photometric redshift 
estimation. Nevertheless, empirical methods do in principle 
allow the best redshift estimation possible for each popula- 
tion. For populations with strong spectral features but lacking 
accurate template spectra, empirical algorithms are expected 
to show marked improvement over template-fitting methods. 

To improve the redshift estimation of z < 0.5 PAH-emitting 
starburst galaxies, as well as very active galaxies, for which 
the AGN component dominates the observed mid-IR emis- 
sion, the AGES spectroscopic sample was used to train a 
neural net algorithm. We adopted the public code ANNz 
(Collister and Lahav 2004) for this purpose. 

Briefly, the ANN is trained to match a set of observational 
inputs (in this case galaxy photometry) to a set of known out- 
puts (the spectroscopic redshifts) by minimizing a cost func- 
tion. The form of the cost function is determined by the ar- 
chitecture of the ANN, which specifies the number of hid- 
den nodes betwee n the input and output nodes. Following 
iFirth et alJ (120031) . we employ a 7:10:10:10:1 architecture, 
which takes 7 inputs, the fiu/fl/[3.6][4.5][5.8][8.0] photom- 
etry, has three 10-node hidden l ayers and a single output 
node, the photometric redshift (see Collister and Lahav 2004; 
IFirth etaljr2003HLahavetalJll996l for details). The inclu- 
sion of the [5.8] and [8.0] photometry (or limits) allow the 
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PAH features in starburst galaxies and the mid-IR excesses 
in active galaxies to be fitted, or at least differentiated. Only 
objects observed in all of the above bands were used in the fit; 
the near-IR data covering only a subset of the field was not 
used. 

The AGES sample was divided into training, validation and 
testing subsets, containing 7500, 2500, and 5052 objects, re- 
spectively. The cost function is fitted on the training set and 
evaluated on the validation set after each iteration. When 
trained, the algorithm can be run on the independent testing 
set to determine the accuracy of the method. The "consen- 
sus" median prediction from a committee of ANNs, each the 
result of independent training sessions (initialized with differ- 
ent random seeds), prod uces the mos t reliable redshifts and 
error estimates (e.g., Firt h et aLll2003l) . Although the current 
implementation of ANNz does not incorporate the photomet- 
ric errors in its determination of the photometric redshift, it 
does use them to estimate the redshift uncertainties. A com- 
mittee of 10 ANNs was adopted in this work. 

The results are plotted in Figure |2] for both the full AGES 
sample (top panel), and the smaller test set which was inde- 
pendent of the ANN training. The results for these two sam- 
ples are essentially identical, indicating that our training sam- 
ple is large enough to span the distribution of spectral types 
and redshifts in the AGES sample. The dispersion over all 
redshifts, a ~ 0.13, is dominated by the z > 1 AGN, whose 
weak broad-band spectral features lead to greater redshift un- 
certainty for any algorithm. Considering just the z < 1 AGES 
galaxies, we find the dispersion drops to a ~ 0.08, whereas 
for the z > 1 AGN it rises to a ~ 0.4. This is typical of the 
accuracy of oth er attempts to measure AGN redshifts photo- 
metrically (e.g., Babbed geetai]|2004l iKitsionas et alJl2005l 
IWeinstein et al.l2004l). 
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FIG. 2. — ANNz photometric vs. spectroscopic redshift for the full AGES 
sample (top panel), and the testing sample alone (bottom). The clipped re- 
sults are quite similar in both samples, with very accurate (<r ~ 0.05) redshift 
prediction for the normal galaxy sample at z < 1 (including PAH-emitting 
galaxies), and more modest accuracy (a ~ 0.3—0.4) at z > 1 where the sam- 
ple is dominated by QSOs and AGN. 



3.4. Hybrid Approach 

In assessing the advantages and drawbacks of the two red- 
shift algorithms presented so far, an obvious complementar- 
ity is apparent. The ANN provides excellent redshift esti- 
mation for low-z galaxies bright enough to have well mea- 



sured photometry in all 4 IRAC bands, including strongly 
PAH-emitting galaxies. Furthermore, the ANN redshift ac- 
curacy for AGN is superior at all redshifts to that achievable 

by template-fitting methods. Due to the IR-excess of these 

i i i I 

active galaxies (Stern et al. 2005), large numbers of them are 

well-measured in all four IRAC bands, at least out to z ~ 2. 
Therefore, within the bright limits of the AGES survey, these 
two populations are well represented and hence well cali- 
brated in the ANN. 

The template-fitting method, while not reliable for 
accretion-dominated active galaxies, is quite robust for nor- 
mal galaxies in the z < 1 AGES sample, barring low redshift 
strong PAH-emitters. This method has the critical advantage 
that it is expected to be robust outside the parameter space 
in which it was tested. In particular it should be reliable to 
magnitudes much fainter that the AGES sample, and out to 
higher redshift. Another important advantage of this method 
is that it allows straightforward generation of redshift prob- 
ability functions, which are essential for many applications. 
Finally, it produces restframe properties for the fitted galax- 
ies, such as absolute magnitudes. 

We therefore construct a hybrid photometric redshift 
sample, combining the above methods according to their 
strengths. While it is tempting to simply adopt the ANN 
method for all bright galaxies, we limit its sphere of influ- 
ence to the color-selected AGN and PAH-emitters as defined 
below. 

Selection of objects for ANN redshift estimation was made 
using color cuts in [3.6] -[4.5] vs. [5.8] -[8.0] color-color 
space, illustrated in Figure [3] Model tracks illustrate where 
representative quiescent and starburst galaxies appear from 
< z < 2 in this plot . Obje cts lying in the "AGN wedge", 
defined in iStern et aTl d2005l) . were classified as likely ac- 
tive galaxies. On the other hand, objects outside the AGN 
wedge and with [5.8] -[8.0] > 1 were taken to be potential 
PAH galaxies. For both populations we only estimated red- 
shifts for objects matching the IRAC flux limits of the full 
spectroscopic sample, which were [17.5, 17.0, 16.2, 15.5] for 
AGN and [16.75,16.6,16.2,14.25] for PAH-emitters in the 
3.6, 4.5, 5.8, and 8.0 micron bands, respectively. The limits 
are generally deeper for the active galaxies since unresolved 
AGN candidates were targeted to a greater depth in AGES. 
No extrapolation to fainter magnitudes was allowed for either 
object class. In addition, adequate photometric coverage was 
required in all 7 bands that were used in the calibration of the 
ANN. The strong flux constraints, in particular in the [5.8] 
and [8.0] bands, limit the sample to the very brightest AGN 
at z < 2.5 and PAH-emitters at z < 0.5, as in the AGES sam- 
ple. The final selection included 3681 AGN and 4766 PAH- 
emitters for which the template-fitting redshifts were replaced 
with ANN redshifts. Representing only ~ 4% of the galaxy 
sample, this approach serves primarily to reduce the number 
of outliers. 

3.5. Comparison to Spectroscopy 

In this Section we demonstrate that the hybrid method sur- 
passes the simple template-fitting algorithm in redshift accu- 
racy for both galaxies and AGN, though the improvement is 
more substantial for the latter. We adopt the QSO/AGN tar- 
geting criteria employed in the AGES survey to distinguish 
between normal and active galaxies. These criteria, described 
fully in Kochanek et al. (in preparation), identify AGN by 
combining optical morphology with flux or color cuts in x- 
ray through radio wavelengths. Figure0]compares the hybrid 
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FIG. 3. — IRAC color-color plot showing the color criteria for inclusion 
in the neural net component of the hybrid photometric redshift algorithm. 
Grey contours illustrate the colors of the full photometric redshift sample. 
The likely AGN and PAH-emitters are overplotted in red and blue points, 
respectively. In addition to these color selections, objects are subject to strict 
AGES-defmed flux limits, as defined in the text. Mo del tracks of two < 
z < 2 non-evolving galaxy templates from Devriendt et al. 1 1999) are also 
plotted. The solid curve represents M82, a starburst galaxy; the dashed curve 
shows NGC 4429, an SO/Sa galaxy with a star formation rate approximately 
4000 times lower. The tracks begin at z = at the diamond symbols, and 
circular annotations are made at z = 0.25, 0.5, 0.75, 1.0, 1.5, and 2.0. 



(top panel) and template-fitting (bottom panel) photometric 
redshift predictions to the spectroscopy for the AGES sam- 
ple of galaxies (dots). The difference between the methods is 
subtle, but most evident at low redshift (z ~ 0.2) where the 
improvement for starburst PAH-emitters is apparent. Also 
plotted is the in-house spectroscopic galaxy sample (filled 
squares), which extends to much fainter flux limits (R ~ 25) 
and higher redshifts than the AGES galaxies. As such it pro- 
vides a valuable, independent test of the method. Clearly red- 
shift estimation is reliable to z ~ 1.5 for this deeper sample. 
There are a handful of outliers between 0.5 < z <j 1, and evi- 
dence of possible systematic errors at z < 0.25 where the Bw 
filter is not fully blueward of the 4000 A break. 

Various measures of the accuracy of the hybrid and 
template-fitting methods are given in Tabled The basic result 
is that a redshift accuracy of a « 0.06 (1 +z) is being achieved 
for 99.5% of galaxies in the AGES sample. For the more chal- 
lenging in-house sample, 98.5% of galaxies have a redshift 
accuracy of a w 0.10(1 +z). The hybrid method also reduces 
the number of outliers, as seen in the fraction of 3<r outliers. 
The final two columns of Table ^report the dispersions for 
subsamples constrained to contain 95% of the objects, thereby 
allowing a direct comparison of the hybrid and template meth- 
ods. Clearly the hybrid method is superior, reducing the dis- 
persion by over 20% compared with standard template-fitting. 

Figure |3 shows a similar comparison for the AGN sam- 
ple. Here the improvement of the hybrid technique over sim- 
ple template-fitting is dramatic, though not surprising. With 



no AGN template in the mix, the template-fitting algorithm 
should not be expected to succeed for accretion-dominated 
objects. This is clear from the lower panel of this figure, 
where beyond z > 1 the results are essentially useless. The 
fairly accurate results at low redshift are likely for galaxies 
which, though active, have luminosities dominated by fusion 
processes. In the bright AGES sample, these galaxies would 
only be visible at modest redshifts, whereas the z > 1 sample 
should be almost entirely composed of extremely luminous 
quasars (some of which are also present at low redshift). 

As a check on this hypothesis, histograms of the mor- 
phological stellarity indicator discussed above are plotted in 
Figure [6] for two AGN subsets, split according to their ob- 
served redshift accuracy and range. Those AGN for which 
the template-fitting algorithm works reasonably well, taken 
to be those with both spectroscopic redshifts and absolute red- 
shift differences of less than unity, are plotted as the solid his- 
togram in Figure[6] The complement, for which the algorithm 
largely fails, is plotted as a dot-dashed histogram. Thus the 
bimodality observed in the lower panel of Figure|5]is strongly 
reproduced in the morphological stellarity measurements in 
Figure [6] The AGN for which the template-fitting method 
works are clearly resolved objects for which the nuclear emis- 
sion does not dominate the flux of the galaxy. Conversely, the 
redshift failures are overwhelmingly unresolved sources, con- 
sistent with the expectation of nuclear-dominated emission 
from very luminous AGN and QSOs. In marked contrast with 
the template-fitting redshifts, the hybrid redshifts (top panel of 
Figure[5} are quite accurate for all AGN to z ~ 2, with no obvi- 
ous systematic issues or significant occurrence of catastrophic 
errors. Beyond z > 2 the redshifts for AGN are systematically 
underestimated, presumably due the relative paucity of cali- 
brators at these high redshifts. 

Statistics of the redshift accuracy for AGN are given in Ta- 
ble |2 The hybrid redshift dispersion is<r«0.14(l+z) for 
over 97% of the active galaxies. The tremendous improve- 
ment achieved by this method, apparent in Figure|5] is borne 
out in a direct statistical comparison of the 95%-clipped sam- 
ples. The dispersion of the hybrid method is a factor of 3 
smaller than that achieved with galaxy template fitting. 

3.6. Comparison with an Independent Photometric Redshift 

Catalog 

We have also verified that the hybrid photometric redshift 
algorithm presented here is in excellent agreement with an 
independent photometric redshift catalog in the same field. 
Brown et al. (in preparation), adopting a pure neural net ap- 
proach, have generated a photometric redshift catalog from 
independently extracted multi-color catalogs, photometered 
with an original code. By making fainter copies of the AGES 
spectroscopic galaxies, they have effectively extended the cal- 
ibration set to much fainter magnitudes. In addition to the 
photometric data, structural information, in the form of sizes 
of the major and minor axes, was also incorporated into the 
neural net for bright objects to improve accuracy at low red- 
shift. Using an independent calibration of the neural net on 
this extended sample, they derive neural net photometric red- 
shifts for a large optically-selected sample in Bootes. 

At relatively bright magnitudes (R < 23) objects with SEx- 
tractor CLASS_STAR parameters greater than 0.85, mea- 
sured in the best-seeing optical band, are taken be stars (Jan- 
nuzi et al. in preparation) and are removed for this compar- 
ison. The stellar contamination is negligible faintward of 
this limit and no attempt is made to remove fainter stars. 
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FIG. 4. — Hybrid (top) and template-fitting (bottom) photometric redshifts for the AGES galaxy sample (dots) and the deeper, in-house spectroscopic sample 
(filled squares). The hybrid method produces a major improvement in z < 0.5 starburst galaxies, and, as quantified in Tablelll reduces the overall dispersion by 
over 20% compared with the standard template-fitting method. 



Table 1. Photometric Redshift Accuracy for Galaxies 



Sample 


Algorithm 


Undipped 


3ct 


Clipped 




95% Clipped 


a 


<t/(1+z) 


% Rejected 


a 


<t/(1+z) 


a 


T/0-+Z) 


AGES 


Hybrid 


0.143 


0.105 


0.52 


0.077 


0.060 


0.060 


0.047 


In-House 


Hybrid 


0.397 


0.185 


1.46 


0.160 


0.096 


0.101 


0.059 


AGES + In-House 


Hybrid 


0.160 


0.109 


0.53 


0.080 


0.061 


0.062 


0.048 


AGES 


Template 


0.230 


0.170 


0.60 


0.102 


0.081 


0.079 


0.061 


In-House 


Template 


0.498 


0.253 


2.51 


0.190 


0.111 


0.127 


0.081 


AGES + In-House 


Template 


0.245 


0.174 


0.69 


0.104 


0.082 


0.081 


0.062 



Table 2. Photometric Redshift Accuracy for AGN 







Undipped 


3a 


Clipped 




95% Clipped 


Sample 


Algorithm 


0- <j/(1+z) 


% Rejected 


(j 


<r/(l+z) 


a <t/(1+z) 


AGES 
AGES 


Hybrid 
Template 


0.473 0.219 
0.998 0.540 


2.94 
0.60 


0.309 
0.918 


0.138 
0.462 


0.255 0.120 
0.797 0.341 
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FIG. 5. — Hybrid (top panel) and template-fitting (bottom panel) photometric redshifts for the AGN sample. The improvement of the hybrid method is clear 
at all redshifts, dramatically so at z > 1. A statistical comparison of the two methods, given in Table|2| demonstrates that the hybrid method reduces the redshift 
dispersion by 65% over the template-fitting method. 
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Table 3. Dependence on Galaxy Magnitude and SED 
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FIG. 6. — Histogram of the SExtractor CLASS_STAR parameter in the 
best-seeingoptical band of two subsets of AGN isolated from the lower panel 
of Figure [5] This parameter varies continuously between zero and unity, 
where the extreme values correspond to sources which are completely re- 
solved and unresolved, respectively. The solid histogram represents those 
objects for which the template-fitting algorithm is reasonably robust, taken 
to be those objects with |Phot-z — Spec-z < 1 and Spec-z < 1. The dot- 
dashed histogram is for those objects excluded by this selection. Clearly the 
template-fitting algorithm works well for resolved AGN and fails for truly 
quasi-stellar objects. 



Over the redshift range where spectroscopic calibrators ex- 
ist, < z < 1.5, the inter-catalog 95% clipped redshift dis- 
persion for the sample of galaxies common to both samples 
is a = 0.09, or a = 0.05(1 +z). These two sets of redshifts 
were computed using independent galaxy photometry and er- 
ror estimates, and for the vast majority of objects, different 
photometric redshift algorithms (template-fitting vs. artificial 
neural net). Yet there is excellent agreement, similar in ac- 
curacy to that demonstrated versus spectroscopy in the pre- 
ceding section, probing right down to the 13.3/iJy limit. This 
provides strong evidence that both photometric redshift cata- 
logs are free from substantial systematic errors. 

3.7. Dependence on Magnitude and SED 

Comparing the results for the AGES and in-house galaxy 
samples in Table ^ it is clear that the redshift precision is 
lower for the fainter in-house sample. In general the pho- 
tometric redshift precision depends on both photometric S/N 
and galaxy spectral type, with slightly smaller uncertainties 
typically achieved for redder galaxies due to their larger con- 
tinuum breaks. Table [3] quantifies the hybrid redshift accu- 
racy in differential magnitude bins for galaxies classified in 
the template fits as earlier or later than an unevolved CWW 
Sbc galaxy. The available spectroscopy beyond the AGES 
4.5yum statistical limit of 15.7 mag is both incomplete and 
inhomogeneous. The tabulated redshift accuracies therefore 
may not be representative of those which would be obtained 
for magnitude limited samples to fainter depths. 



4.5/jm Mag Range 


Early-Type 


Late-Type 


a 


<r/(l+z) 


a 


<x/(l+z) 


14.5-15.0 


0.039 


0.031 


0.038 


0.030 


15.0-15.5 


0.038 


0.029 


0.046 


0.034 


15.5-16.0 


0.040 


0.029 


0.049 


0.036 


16.0-16.5 


0.050 


0.037 


0.071 


0.054 


16.5-17.0 


0.069 


0.054 


0.083 


0.063 


17.0-17.5 


0.074 


0.060 


0.073 


0.057 



NOTE. — The AGES + in-house galaxy sample is divided 
into early- and late-type samples according to the best-fit 
SED templates; objects with best-fit templates earlier than an 
unevolved CWW Sbc template are classified as early-type, 
whereas those of type Sbc or later are taken to be late-type. 
Statistics are for the 95% clipped sample. 



4. REDSHIFT PROBABILITY FUNCTIONS 

Redshift likelihood functions were constructed by project- 
ing the likelihood surface of redshift and spectral type onto 
the redshift axis. Convolution of these likelihood func- 
tions with a variable-width Gaussian kernel, which increases 
the error with redshift, accounts for the templat e-mismatch 
variance inherent in the method (e.g.. iBrodwin et aTll2006t 
Fernandez-Soto et al. 2002). From the results of Table[fl the 
kernel is taken to be a(z) = 0.06(1 +z). 

This renders the area under these functions reasonable 
proxies of redshift probability density, resulting in approxi- 
mate redshift probability distribution functions (PDFs) for all 
objects in the survey. Due to the excellent redshift accuracy 
for < z < 1.5, no redshift prior was applied to these likeli- 
hood functions, other than the weak prior imposed by the lim- 
ited fitting range. Stronger priors can, of course, be imposed 
for certain kinds of analyses; we employ this methodology in 
computing the galaxy redshift distribution in the next section. 
The statistical validity of the PDFs in the present paper can 
be explicitly confirmed for those objects which have spectro- 
scopic redshifts as follows. 

Redshift confidence intervals are derived from the PDFs by 
associating area with probability density as illustrated in Fig- 
ure We define the la, 2a, and 3a confidence intervals as 
those redshift regions which enclose the top 68.3%, 95.4% 
and 99.7% of the normalized area under the PDFs. While the 
full PDF should in general be used for statistical analyses of 
the galaxy population, confidence intervals defined in this way 
offer a straightforward test of the PDFs. In Table|4]we report 
the fraction of objects which agree with the spectroscopic red- 
shift at each confidence level, for both the full spectroscopic 
sample and for the galaxy sample alone (i.e., excluding ob- 
jects identified as QSOs and AGN in AGES, but including 
PAH-emitting starburst galaxies). We only consider objects 
here with photometry of sufficient quality to be included in 
the m ain photometric redshift sample as described in Section 
l2~2~2l 

The 1 a confidence intervals are approximately Gaussian for 
the full sample in that the spectroscopic redshift is included 
in this photometric redshift confidence interval for 70% of the 
objects. The 2a and 3a confidence intervals have inclusion 
rates slightly below the Gaussian expectation, by 3.9% and 
2.7%, respectively. This small outlier fraction is reduced to 
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FIG. 7. — Redshift PDFs for 3 sample galaxies, with the lcr (68.3%), 2cr (95.4%), and 3cr (99.7%) confidence intervals denoted by horizontal lines under the 
curves. 



0.8% and 1.2% when QSOs and AGN are excluded, though 
the lcr interval is slightly conservative in this case. That 
the inclusion of a large fraction of active galaxies affects the 
statistics so little is an indication of the robustness of the red- 
shift probability functions. The lack of strong continuum 
spectral features results in quite broad redshift probability 
functions for these objects, reflective of the larger uncertainty 
in their redshifts. These lead to wide, but valid, confidence 
intervals, which can be incorporated in large statistical stud- 
ies. On the other hand, previous IRA C Shallow survey papers 
jEisenhardt et alJl200H IStern et al.ll2005h have demonstrated 
how these objects can be reliably removed, if desired, using 
photometric information prior to redshift fitting. 

5. SCIENCE APPLICATIONS 

A shared primary goal of the IRAC Shallow, NDWFS and 
FLAMEX surveys is to study structure formation and evolu- 
tion at 0.5 < z < 2, both in the field and in cluster environ- 
ments. In this Section we present two specific large scale 
structure science applications enabled by the hybrid photo- 
metric redshifts and redshift probability distributions derived 
in this paper: a measurement of the 4.5 /im galaxy redshift dis- 
tribution, and the discovery of a high-redshift (z > 1) galaxy 
cluster. 

5.1. Redshift Distribution at4.5[wi 

A key issue in structure formation mo dels is the mass as - 
sembly history of massive galaxies ( e.g.. iFaber et al.] l2006). 
Recent work (Bunker et al. 2006; Mobasher et alJ 120051 
lYanetal.ll 2005) has indicated that massive galaxies form the 
bulk of their stars soon after reionization, and evolve pas- 
sively over most of the h istory of th e universe (tBundv et al.l 
120051 lJuneau^TaIll2005l iTreu et al.ll2005l) . While hierarchi- 
cal models can accommodate a modest number of early, mas- 
sive halos, these results were not predicted in advance of the 
observations. 

More generally, any viable galaxy formation theory must 
predict the low order moments of the mass distribution, in- 
cluding the redshift distribution and the autocorrelation func- 
tion. Previous attempts to constrain galaxy formation mod- 
els using these moments have been limited by the difficulty 
in rel ating optical light to mass, the natural theoretical vari- 
able. [Kauffman n and Chariot! ( 119981) suggested employing a 
K-hwA selection to minimize this source of uncertainty, and 



ICimatti et alJ (120021 hereafter K20) subsequently presented 
the K s < 20 redshift distribution. Their sample had a median 
redshift of z mec j ~ 0.80, and a distribution that was best re- 
produced by pure luminosity evolution models. 

However, the /T-band is only a good proxy for stellar mass 
at redshifts where it samples the rest-frame near-IR stellar 
peak. At 1 < z < 2, where significant elliptical galaxy as- 
semb ly is expected in hierarchical models (e.g., Kauffmann 
1996), the stellar peak is firmly in the Spitzer/IRAC bands. 
We therefore present in Figure [8] the 4.5/im galaxy red- 
shift distribution derived from the 13.3/iJy sample described 
above. The survey is over 85% complete to this limit, based 
on the recovery fraction of artificial stars in standard com- 
pleteness simulations. Very strict multi-band masking was 
employed in determining the redshift distribution, rejecting 
areas not containing valid coverage in all of the key optical 
(B W RI) and IRAC ([3.6][4.5]) bands. This resulted in a fi- 
nal unmasked area of 7.25 deg 2 . Stars are rejecte d at b right 
magnitudes only, using the criteria described in 33.61 The 
redshift distribution is calculated in two ways. The simple 
histogram is derived using the best-estimate redshifts of the 
hybrid method. The curve in the figure shows the result of 
summing up the full normalized redshift probability func- 
tion fo r each galaxy, following the method of Brodwin et al. 
(2006). The prior in this method is taken to be the approx- 
imate N(z) estimated from a direct summation of the galaxy 
likelihood functions. The vertical dashed line represents the 
redshift limit of the available spectroscopy; beyond this limit 
the robustness of the photometric redshifts has not been ex- 
plicitly demonstrated. The median redshift is z met j = 0.98 for 
the simple hybrid-method histogram, and z me( j = 0.99 for the 
PDF summation method. Assuming the galaxies with for- 
mal photometric redshifts of z > 2 are distributed according 
to the plotted < z < 2 distribution, the median redshifts are 
z med = 0-85 and z mec i = 0.88, respectively. The two methods 
are in excellent agreement, which, while not tautological, is 
nevertheless expected. It provides a measure of confidence in 
the consistency of the methods employed here. 

Despite being shallower than the K20 survey, the IRAC 
Shallow survey has a higher median redshift, owing to a ben- 
eficial, negative K-correction (e.g., Eisenhardt et al. in prepa- 
ration). This deeper reach, along with a 500-fold increase in 
area over the K20 survey, enables stronger constraints to be 
placed on models of galaxy evolution. We will be examining 
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Table 4. Confidence Level Statistics 
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Expectation 
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Confidence Interval 
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< la 


68.3% 


10870/15530 


70.0% 


9722/13043 
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< la 
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12848/13043 
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FIG. 8. — Differential 4.5/im-selected photometric redshift distribution in Bootes. The hybrid redshift catalog was used for the histogram. The smooth curve 
was constructed from summing up the individual galaxy redshift probability functions as described in the text. The dashed vertical line indicates the maximum 
redshift at which the photometric redshift method has been verified against spectroscopy. 



this in detail in a future paper. 

5.2. A High Redshift Galaxy Cluster Search 

We present some early results of a high redshift (z > 1) 
galaxy cluster search underway in the Bootes field. The de- 
tection techniqu e, described in Eisenhardt et al. (in pr epara- 
tion, see also Elsto n et afll2006t IStanford et al.1120051 Gon- 
zalez et al. in preparation), implements a wavelet search algo- 
rithm tuned to identify structure on cluster scales (~ 500 kpc). 
The redshift probability functions are the input to the wavelet 
code, and so the cluster detection is principally dependent on 
the accuracy and statistical reliability of the redshift PDFs. 
The method is independent of the strength, or even presence, 
of the cluster red sequence, and therefore provides an unbi- 
ased window on the era of cluster formation. 

Four cluster candidates were targeted spectroscopically in 



the first half of 2005, and all four were confirmed to be 
z > 1 galaxy clusters, at redshifts ranging from (z) = 1.11 to 
(z) = 1 .41 . The latter c luster, the highest yet f ound in a cluster 
survey, is presented in Stanford et al. (2005). In this Section 
we present one of the newly discovered clusters, a filamentary 
cluster at (z) = 1.24. 

The cluster ISCS J1434.5+3427 is shown in Figure^ The 
cluster-finding power of the Spitzer/IRAC imaging is clearly 
demonstrated in this figure. Though no structure is apparent 
in the optical bands, a striking filamentary structure emerges 
beyond rest-frame 4000 A in the 4.5/im band. The spectro- 
scopic cluster members are circled in red on the greyscale 
4.5/im image. Objects which have the cluster systemic red- 
shift of (z) = 1.24 within their ltr confidence intervals are 
marked by blue squares. 

ISCS J 1434.5+3427 was observed spectroscopically in 
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2005 February with Keck/ LRIS JOke et aljfl995h and 2005 
May with Keck/DEIMOS JFaber et alJ 120031) . For a more 
detail ed descr i ption of these spectroscopic observation see 
IStanford et alJ (120051) and Desai et al. (in preparation). The 
eight spectroscopic members within Az = 0.01 of the sys- 
temic redshift confirm the reality of this cluster at (z) = 1.241. 
Details of these members are given in Table [5] An estimate 
of the cluster velocity dispersion is deferred until additional 
spectroscopy yields more members. In addition, deep follow- 
up imaging observations with HSTI ACS and SpitzerflRAC are 
underway. 

6. SUMMARY 

Accurate photometric redshifts, calibrated using over 
15,000 spectroscopic redshifts, have been computed for a 
4.5^m sample of 194,466 galaxies in the 8.5 deg 2 IRAC Shal- 
low survey. A hybrid technique, in which a standard tem- 
plate fitting code is augmented using a neural net approach, 
was adopted to optimize the redshift accuracy for both ac- 
tive and low-z starburst galaxies without compromising re- 
liability for the general galaxy population. This is primar- 
ily enabled by the fact that these two populations, which are 
often troublesome for template-fitting codes, are particularly 
well represented in the AGES sample. Between < z < 1.5 
the resulting hybrid algorithm has a demonstrated accuracy of 
a < 0.06(1 +z) for 95%, and a < 0.1 (1 +z) for 98.5%, of the 
galaxy population. For over 97% of the active galaxies, the 
redshift accuracy is a < 0.14(1 +z). 

Redshift probability functions have been computed for all 
objects directly from the template-fitting algorithm. Com- 
parison with the large spectroscopic sample has verified the 
statistical validity of these functions, and in particular, the 
reliability of confidence intervals derived from them. These 
confidence intervals, or indeed the full probability functions, 
can be reliably used in statistical studies of the galaxy popu- 
lation. Several such programs are underway, and we present 
in this paper two new results which employ them. 

The 4.5/im-selected galaxy redshift distribution, a primary 
observable for confronting theories of structure formation, 
was computed using both the hybrid photometric redshifts and 
the full redshift probability functions. The methods yield en- 
tirely consistent results. This measurement is provided in an- 
ticipation of future model predictions extending into the mid- 
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IR, where flux is closely related to stellar mass above z > 1. 

Another program making extensive use of these redshift 
PDFs is a search for high redshift (z > 1) galaxy clusters. 
We presented one such cluster, ISCS J 1434.5+3427, spec- 
troscopically confirmed at (z) = 1.24, which was discovered 
by incorporating the redshift PDFs in a wavelet search algo- 
rithm. Spectroscopic confirmations of two similarly discov- 
ered high redshift clusters, a t (z) = 1.11 and (z) = 1.41, are pre- 
sented in companion papers ( Elston et alJ2006tlStanford et alJ 
120051) . A complete description of the cluster survey sam- 
ple and methodology is presented in Eisenhardt et al. (in 
preparation), along with the spectroscopic confirmation of a 
(z) = 1-37 cluster. 
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FIG. 9. — 5' x5' Bw, I, [4.5], and composite color image of cluster ISCS J1434. 5+3427 at (z) = 1.24. North is up and East is left. Notice the striking filamentary 
morphology that emerges in the 4.5/im band, and is quite prominent in the composite image. Spectroscopically confirmed cluster members are denoted by red 
circles on the [4.5] image. Objects which contain the cluster systemic redshift within their lcr photometric redshift confidence levels are denoted with blue 
squares. 



Table 5. Summary of Spectroscopic Members 













95% Confidence 










ID 


R.A. a 


Dec." 


[4.5] b 


phot-z 


Interval 


spec-z 


Sz 


Date 


Instrument 


IRAC J143430.3+342712 


14:34:30.36 


+34:27:12.1 


15.81 


1.26 


[0.98, 1.63] 


1.2365 


0.0005 


UT 2005 Feb 10 


LRIS 


IRAC J143428.6+342557 


14:34:28.66 


+34:25:57.7 


15.19 


1.12 


[0.89, 1.43] 


1.238 


0.003 


UT 2005 Feb 10 


LRIS 


IRAC J143421. 6+342656 


14:34:21.64 


+34:26:56.2 


16.67 


1.09 


[0.84, 1.37] 


1.2502 


0.0005 


UT 2005 Feb 11 


LRIS 


IRAC J 143429.2+342739 


14:34:29.21 


+34:27:39.1 


16.99 


1.06 


[0.68, 1.49] 


1.2436 


0.0005 


UT 2005 Feb 11 


LRIS 


IRAC J143430. 1+342657 


14:34:30.15 


+34:26:57.2 


15.90 


1.17 


[0.92, 1.48] 


1.23 


0.01 


UT 2005 Feb 11 


LRIS 


IRAC J143428.0+342535 


14:34:28.05 


+34:25:35.8 


16.64 


1.19 


[0.89, 1.57] 


1.240 c 


0.002 


UT 2005 May 07 


DEIMOS 


IRAC J143418.4+342733 


14:34:18.41 


+34:27:33.1 


16.42 


1.21 


[0.96, 1.52] 


1.240 


0.001 


UT 2005 May 07 


DEIMOS 


IRAC J143430.6+342757 d 


14:34:30.63 


+34:27:57.2 


15.35 


0.95 e 


[0.00, 4.17] c 


1.242 


0.002 


UT 2005 May 07 


DEIMOS 



"Coordinates are J2000. 

b Vega magnitude at 4.5/xm; mag = 179.5 Jy. 

c Lower quality redshift due to sky lines superimposed on [O II] feature. 
d MIPS source. 

e Photometric redshift for this object is from the neural network. 
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